Optimal Polynomial Admissible Meshes on Compact 
Subset of with Mild Boundary Regularity^ 



It has been proved [6] that discrete least squares polynomial approxi- 
mation performed on Polynomial Admissible Meshes, say AM, enjoys nice 
property of convergence. Optimal AMs [H] are AMs which cardinality 
grows with optimal rate w.r.t. the degree of approximation. Here we present 
two new results in this framework. 

In Section 2 we show that any compact subset of M*^ that is the closure 
of a bounded star-like domain preserving a uniform interior ball condition, 
say UIBC (cfr. [T]), admits an optimal AM. This extends a result of A. Kroo 
[H] proved for star shaped domains and is closely related to the recent 
preprint [T5] . 

In Section 3 we prove constructively the existence of an optimal AM in any 
^^'^ compact subset of M'^, this is done recovering a multivariate counterpart 
of celebrated Bernstein Inequality via the distance function. 

Keywords: Admissible Meshes, Multivariate Polynomial Approximation, 
Uniform Interior Ball Condition, ^^'^ domains. Distance Function. 



1. Optimal Polynomial Meshes 

Let K be any compact polynomial determining subset of M.'^ and {v4„}n 
a sequence of subsets of it. An is said to be a n-degree polynomial admissible 
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mesh (AM) of constant (i.e. not depending on n) C > if CardA„ = 0(n*) 
for a suitable s G M and if for any polynomial p of degree at most n in o? 
variables {p G ^"(M'^)) the following inequality holds true 

IHk<^^IHlA„. (1) 

Here and throughout the paper we let := sup^.^^^ 1/(3^)1- 

If instead we allow C to grow at most polynomially w.r.t. n (i.e. there 
exists s' > such that C„ = (9(n* )), then An is said to be a weakly admissible 
mesh (WAM)[n]. 

The most relevant motivation to introduce such definitions is that, under 
mild regularity assumptions on a target function /, one has the uniform con- 
vergence on K of discrete least squares polynomial approximation performed 
samphng f on An- More precisely in [H] Authors proved that 

11/ - AaJWk <{l + Cn (||/|k(l + ^/CaTd{An)))) dn{f,K), (2) 

where Aa_„ is the discrete least squares polynomial projector, Cn is the con- 
stant of the n degree WAM An and dn{f,K) := infpg^n(]Rd) ||/ — 

An Admissible Mesh can be built [51 Theo.5] by a grid thick enough on 
any compact set satisfying a Markov Inequality for polynomials [5], that is 
there exists a (Markov) constant M > and a (Markov) exponent r > 1 
such that for any polynomial p G ^"(M'^) we have 

WI'^pIooW K <Mn'-\\p\\K. (3) 

Here | ■ |oo is the maximum norm in M"'. The existence of a Markov Inequality 
holding on a given fat (i.e. intRd K = K) compact set K is not too restrictive. 
For instance in the case of real variables one has a Markov Inequality with 
exponent 2 for any compact set satisfying a Uniform Cone Condition^ thus for 
the closure of any bounded Lipschitz domain. However ^ is preserved with 
an exponent greater than 2 even by sets admitting also cusps of polynomial 
type, namely Uniformly Polynomially Cuspidal (UPC) sets, we refer to [20] 
and [21j for all involved definitions and a specific treatment on this topic. 

WAMs and AMs were shown to enjoy nice properties as stability un- 
der unions and tensor products, stability under affine mappings, WAM can 
be easy computed on sets which are the polynomial images of a set where a 
WAM is known, moreover any set of good interpolation points forms a WAM. 
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Hereafter we refer to any unisolvent set F„ of interpolation points for poly- 
nomials of degree n having slowly increasing (i.e. limsup„(Le6(F„))^/" = 1) 
Lebesgue constant as a good interpolation set. 

For a survey on WAMs we refer to [16j. After this paper we proved two 
interesting results concerning mapping and perturbations of WAMs |22] |21] ■ 

A striking property of WAMs is that a specific good interpolation set, 
say approximate Fekete points (AFP), see for instance [TB], can be extracted 
from a WAM by standard numerical Linear Algebra. These sets of points 
are an approximate solution to the problem of finding true Fekete points 
(i.e. := dim^"'(M'^) = ("^'^) points that maximize the Vandermonde 
determinant) of the given compact K in the following sense. The uniform 
probability measures of true and approximate Fekete Points have the same 
weak* limit [21 Theo. 1], the Pluripotential Equilibrium Measure [25], see [TT] 
for a remarkable and more accessible survey. 

Such an approximate solution is computed by the celebrated QR factor- 
ization with pivoting, the complexity of this procedure grows polynomially 
with the size of the problem, this is one of the reasons because holding the 
cardinality of the starting mesh is important. [3] is an example of a successful 
approach in this sense. 

From the inequality (|2| it is evident that a good discrete least square 
approximation scheme performed on a WAM must ideally have slowly in- 
creasing (or even constant) C„ and slowly increasing CardA„. 

On the other hand from the definition of WAM it follows that a WAM 
is an unisolvent set of degree n. Therefore one immediately has the lower 
bound 



Thus it is rather natural to introduce the definition of optimal admissible 
mesh as any AM having optimal (i.e. 0{n'^)) cardinality [14J. 

In |jl4J the Author refers to any star-shaped compact set K as if its 
Minkowski functional has a Lipschitzian gradient, i.e. 



for a suitable positive C. Under such assumption he shows that K possesses 
an admissible mesh Yn of degree n with 



Card An > dim ^ 



n 




V^k{x) - V^K{y)\ <C\x-y 



a 



CardF„ = 0{n-^ 



(4) 
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In particular he notices that this assures the existence of optimal AMs for 
convex bodies and the closure of any star-shaped bounded domain. 
In the meanwhile of writing this paper we received a new preprint by A. 
Kroo where ^5] Theorem 3] the Author improves his estimate |4] by a smart 
use of Minkowski functional smoothness, he shows that such a new estimate 
is also sharp. 

In [ll] he also conjectured that any real convex body has an optimal 
admissible mesh. 

In this work we build such optimal admissible meshes on two relevant 
classes of compact sets. 

In section 2 we work on star shaped compact sets in M'^ with as less as we 
can boundary regularity assumptions, namely we consider the closure of any 
bounded domain Q such that a Uniform Interior Ball Condition (UIBC) is 
enjoyed by fl, thus we extend the result proved in [23] to the case d > 2. 

For the reader convenience we recall that a set f2 G M.'^ is said to enjoy a 
UIBC of parameter (or radius) r > iff for any x G dQ there exists u G S"'^^ 
such that B{x + ru, r) C Q. 

In section 3 we deal with compact sets being the closure of bounded 
connected ^^'^ domain. We prove constructively the existence of AMs on 
such class of sets. 

Clearly the two classes have non-trivial intersection, however we prefer to 
give separate proofs under different assumptions because we developed two 
techniques: in the first case we use the star-property of the compact, while 
in the second the distance function. 



2. Optimal AM for Star-Shaped Sets Satisfying a UIBC 

In Approximation Theory it is customary to consider as mesh parameter 
the fill distance h{Y) of a given finite set of points Y w.r.t. a compact subset 
X of R"^. 

h(Y) := sup inf \x — y\. (5) 

In this definition it is not important whether the segment [x,y] lyes in X 
or not. If one wants to control the minimum length of paths joining x to y 
and supported in X then he should consider the following straightforward 
extension of the concept of fill distance above. 
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Definition 2.1 (Geodesic Fill-Distance). Let y be a finite subset of tlie 
locally complete set X C M^, then we set 

i44,y(X) := {7 : supp(7) C X,7(0) = x,7(l) = Var[7] < 00} 

and define 

hx(Y) := sup inf inf Var[7], (6) 
the geodesic fill distance of Y over X. 

Here and throughout the paper we denote by Var[7] the total variation of 
the curve 7, 

N 

Var[7] := sup sup \'j{ti) - 7(ti-i)|. 

ATgN 0=to<ti-<tN=l^ 

Notice that the locally completeness of X ensures the existence of a length 
minimizer in the class ^x,y '■= {l '■ 7(0) = a;, 7(1) = Var[7] < 00} provided 
it is not empty, that is if there exist a rectifiable curve ip connecting any x 
and y in X such that Var['?/'] < L < 00, thus if X has finite geodesic diameter 
we can replace inf^g^^ Var[7] by min^g^^^ Var[7] in 

Now we want to build a mesh on the boundary of a connected bounded 
domain satisfying a UIBC prescribing the geodesic fill distance and control- 
ling the cardinality of the mesh, then we use such a "geodesic" mesh to build 
an optimal AM for the closure of the domain. 

We base our construction on a modification of the following proposition 
proved in [12], for a complete treatment on rectifiablity we mention the out- 
standing monograph [TT] . 

In [12], Prop. 2.4] authors deal with cone condition at the boundary of 
parameters r,a, that is Vx G dQ 3i> := z/(x) e S'^"^ such that the open 
(positive half of the) cone of width a and radius r having the vertex in x lies 
in Q: 

X + K{a, r, z/) C , where 

K{a,r,iy) := {y E : y ■ > cosa}. 

Proposition 2.1. [12] Let K dW^ he a compact set satisfying an interior 
cone condition at the boundary of given parameters r,a > 0, then dK is 
^rf-i rectifiable, and is contained in a finite union of Lipschitz graphs. 
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We stress that the following is an easy corollary that follows from the 
previous proposition. Anyway we prefer to give a detailed proof for the sake 
of completeness and to be clearer in the construction. 

Proposition 2.2. Let fl be a bounded set in M"^ satisfying the interior cone 
condition at the boundary with parameters a,r > 0, then for any h > there 
exists {Yn}n G X := dfl 

(i) cardy„ = e»((f)'-') 

(ii) hx{Yn) < h/n. 

Proof. Notice that X := dfl is closed and bounded, hence compact, more- 
over 

dfl = lyj^esrf-i {{x :x + rK{a/2, r, u) (zU}n dfl) , 

thus by compactness there exist a natural number M and a finite sub-covering 
of cardinality M: 

dQ - y,e{i,2,...,M} {{x:x + rK{a/2, r, z/,) C H} n dn) := ^e{i,2,...,M}^j 

Let us now show that any Xj is a finite union of Lipschitz graphs, for let us 
fix j and rotate coordinates by G SOa such that Rii'j) ■ I'j = Cn, more- 

over take a finite covering of Xj made by coordinates-cubes Q^j,Q'j, . . . ,Qj 
having diameter less or equal to r. 

Let us fix s e {1, 2, . . . , 5"} and x,y e Qj nXj, by the very construction 
we have 

X ^y + {R{vj)K{a, r) U R{-Uj)K{a, r)) 

and the same holds true exchanging the role of x and y. For take any y e 
R{—i'j)K{a,r) and assume by contradiction y e Xj then Xj 3 x & y + 
R{vj)K{a,r) C Q, and this is a nonsense. 

Now we write x = {x', Xn),y = {y', yn) and consider the line segment [x, y\ 
and its {d — l)-projection [x\ y'], by elementary geometry the angle between 
this two segments is bounded in modulus by (tt — a/2)/2 and thus we have 

\xn - yn\ < tan((27r - a)/4)\x' - y'\. 

Therefore Xj fl is a compact piece of the graph of a Lipschitz function 
(fj^s '■ IR'^"^ D HQj — > K, where UQj is the d — 1 coordinate projection onto 
the first d — 1 coordinates of Qj. 
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Let us pick a uniform grid A^i^ C IIQ^ such that /ing^l^^'^) < 



h 



we have 



Let us set F„ := U^l^uf^i Trivially Card(AJ;'") = 0{n'^-^) and hence 
Card(A„) = MSOin"^-^) = 0{n'^-^). 

Now take any x E X then there exist j G {1,2... M} , s G {1, 2, . . . , S*} 
and C, G LIQ^ such that ^Pj,s{0 — ^- Notice that the restriction of a Lipschitz 
map to (the support of) a rectifiable curve is trivially a Lipschitz map having 
at most the same Lipschitz constant. 

Let us consider ^j^g '■= o.rgmin^^j^3,s\^ — rj\ and the segment 7 connecting 
f to E-i , with length at most / , thanks to Area Formula [TTl Theo. 

3.2.3] we have 

Var[7] = ^^(supp7) < ^{l + L^)^=L== < h/n. 

n^{l + L^) 

Therefore 

Now we are ready to state and prove our main result of this section, it 
should be compared to the recent preprint [I5l Theorem 3]. The results 
achieved this new manuscript are more general, still they do not cover the 
case of a domain satisfying UIBC but not being ^^'^'^/rf > 2 globally 
smooth. 

Theorem 2.3. Let Q be a connected bounded star shaped domain in M'^ sat- 
isfying a UIBC, then K := Q has an optimal polynomial admissible mesh. 



Proof. We can suppose wlog the center of the star to be by stabihty of 
AM under euclidean isometries |16j . 

Let us consider a° , a^, . . . , a^", the set of 2r;, + 1 Chebyshev points for 
the standard unit interval and set := r(l + a^) for any r > i = 
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1,2, ...2n + 1. By a well known result ([U]) the set Bnij) of all &^(r)'s 
(varying the index i) is an admissible mesh of degree n and constant \/2 for 
the interval [0, r]: 

lbll[o,r.] < v^lbb„M Vpe^". (7) 

hBn{\x\), 



X 

\x\ 



Let us take any x G X := dK and consider the set Bn{x) 

notice that Bn{x) C K because K is star-shaped. 

One can set Z„ := Ux(^xBn{x) and notice that the restriction of any 
polynomial of degree at most n in c? variables to any segment is an univariate 
polynomial of degree at most n, then thanks to ([t]) we have 



(8) 



Therefore we can reduce ourself to find an admissible polynomial mesh of 
degree n for Z„, in other words we can say that Z„ is a norming set for K. 

By the UIBC satisfied by VL there exist r := r(fi) > and u(x) G S*^"^ 
such that B{x + rz/(x)) C K for any x G dK. 

Let us consider any Lipschitz curve 7 in X, notice that for ^^-a.e. 
X G supp7 Tx'y C TxdB{x + rz/(x)). Since the ball is a compact alge- 
braic manifold a Markov Tangential Inequality of degree 1 holds true on it 
(see [5] and references therein) , moreover the constant of such inequality is 
the inverse of the radius of the ball. 



dp 
dv 



x) 



< -^n\\p\\Bixo,r) yp e ^'^(M") , \/v G T,B{xo,r). (9) 



Let us recall that any Lipschitz curve 7 can be reparametrized by arc-length 
parameter by the inversion of t (—)■ Var[7|[o,t]], obtaining a Lipschitz curve 

7 : [0,Var[7]] ^ X 

Var[7] = Var[7] 
Lip[l] = 1 =,iri-a.e. 1 7' I 

Therefore by the UIBC and g -a.e. in the support of 7 we have 

d{p o 7) 



dt 



it) 



< |Vp(7(t))-7WI< 
7'(t)|n 



< 



n , 



^ \\p\\Bm)+rum)),r) < -\\p\\k- 



(10) 

(11) 
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By Proposition 2.2 we can pick Yn{h) on X such that hx(Yn{h)) < h/n 
and CardYn{h) = 0{n'^~^). Let us pick it with h = ^. 

Let us now pick any x & X and consider 7, the geodesic arc connecting 
the (one of the) closest point yl^ of Yn to x and x itself parametrized in the 
arc- length parameter, we have Var[7] < h/n = 

By Lebesgue Fundamental Theorem of Calculus for any p G ^"(M'^) one 

has 



\p{x)\ < \p{yn)\ + 

< b(l/n)l + 

< \piyn)\ + 



^"■■t^l d{p O 7) 
dp 





Var[7] 
h/n „ 



< 



(0 







Kd^<\piy'n)\ + ^\\P\\K 



where in the last line we used (11). Thus we have 

1„ „ 

\\p\\x < \\p\\y„ + ^WpWk- 

By nice properties of rescaling we have also 



bl,X < \\p\\bl,Y„ + l/2\\p\\b^^K < \\P\W„Y„ + ^WpWk, 



(12) 



for consider the omothety G : M"' — ?■ M'^, where 6(x) : 
inequality (12) for g := p o 6. 
Hence we can state that 



< \\p\\ujZobl.Yr^ + V2|blk =: lblU„ + M 



and write the 



K 



(13) 



Now we can use ([8j) and Z„ D U^"o^n^n to get 

WpWk < V2 (^\p\\x„ + ^WpWk^ 

2V2 



K 



H|x„ = 2(v^+l)|b||x„ 



2-^2 

thus Xn is an admissible polynomial mesh for K. The set X„ is the disjoint 



union of 2n + 1 sets 6^F„,thus 



CardX„ = (2n + l)0{n'^-^) = 0{n'^), 
therefore X„ is an optimal admissible mesh of constant 2(v^ + 1). 
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From an algorithmic point of view an AM built by a straightforward 



application of Theorem |2.3| is rather coarse and should be refined. Informally 
speaking such a collocation technique creates AMs that are clustered near 
the center of the star, wile this seem to have no geometrical neither analytical 
meaning. 

This issue could be partially removed by two modification of the con- 
struction. 

First we can work on each := b^^X and build a geodesic mesh on 
it using a Markov Tangential Inequality of different parameter for each F^. 
Namely, for any arc-length parametrized curve 7 in F^ passing through xq 
we consider a (as big as possible) ball B{xo, 6) such that Tx^OB^Xq, r) ^ Tx^"^ 
and we can notice that we can choose 5 = max{(i(F^, X)/2, 6^r}. Then we 
perform the same construction as in proof above to recover an inequality as 



(12). 



Remark 2.5. In [21] [Corollary 1] AMs are stable under small perturbation 
in Hausdorff distance, more precisely if is an AM of constant C for the 
compact set K satisfying a Markov inequality of parameter M and 6 G 
(0, 6*/d), where 6* solves the equation texp {t/2) — 1 = 0, then any finite set 

An C K such that d^^d(An, An) < Mn'^{l+C) ■~ ^" ^^^^ d_%'d 



is the Hausdorff distance. Therefore the application of Theorem 2^ is well 
promising: if f2 is a Markov compact the computation suggested by the proof 
of the theorem needs not to be performed with high accuracy. For instance 
if the center of the star-shaped set is supposing 

min \x\ = m > 0, 
min X ■ vix) = c > 0, 

where ^{x) is the axes of the interior ball, will suffic^ to prove that K 
preserves a Markov inequality of exponent 2 . Thus in such a case the 
conclusion of [21] [Corollary 1] holds true. 

It has been shown (see [Ij) that ^^'^ domains in M'^ are characterized by 
the so called uniform double sided ball condition, that is, f2 is a '^^'^ domain 



^Let X £ n and Cx '■= x/\x\, for any p e ^"{M.'^) such that \\p\\k < 1 one can bound 
l^^j^l and l^^f^l where ejv G S"^^^ O {cx)'^ by an univariate Bernstein Inequality and 
then derive a Markov muhivariate Inequality by the standard procedure. 
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iff there exists r > such that for any x G dVt there exist v G S*^"^ such that 
we have B[x + rf , r) C i7 and B[x — rv, r) C ZVL. Therefore the following is 
a straightforward corollary of our main result. 

Corollary 2.5.1. LetVL he a hounded star- shaped ^^'^ domain, then its clo- 
sure has an optimal AM. 

It is worth to recall that such domains can also be characterized by the 
behaviour of the oriented distance function of the boundary (i.e. 6n(x) := 
d{x,VL) — d{x,ZVL)). For any such ^^'^ domain there exists a (double sided) 
tubular neighbourhood of the boundary where the oriented distance func- 
tion has the same regularity of the boundary as a manifold, this condition 
characterizes ^^'^ domains too. This framework is widely studied in [Hj and 

m- 

3. Optimal AM for ^^'^ Domains by Distance Function 



As we mentioned above in [M] the Author conjectures that any real com- 
pact sets admits an optimal AM, in this section we prove that this holds true 
at least for any real compact set K which is the the closure of a bounded 
^^'^ connected domain Vt. 

Our strategy goes this way: 

• We work out a Bernstein-like inequality for "normal']^ derivatives of 
polynomials depending on the distance function of the complement of 



the considered domain f2. This is done in Proposition by integrat- 
ing the true (one dimensional) Bernstein Inequality of degree n along 
segments given by the metric projection onto the complement of Vl. 
It turns out that this integrals well define a sequence of functions F„ 
which extend continuously to the closure of the domain. 

We split the domain in two parts, namely a tubular neighborhood VL^ 
( in the relative topology of Vt) of the boundary and its complement in 
Vt say Ks- 



•^Here the word normal must be intended in a metric sense w.r.t. dVl and will be clarified 
later. 
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We build for the first one a norming set w.r.t. the domain (by Proposi- 
tion 3.2) by taking the pre-image of a certain number of equally spaced 
points in the interval [0,max^F„]. 

In the second one we prove a Bernstein-like inequality (see Corollary 



22 ) for the directional derivatives of polynomials. 



We give an upper bound for the tangential derivatives of polynomials 



(Proposition 3.3) when restricted to the above level sets defined by 
pre-imaging. 

We build geodesic meshes on level-sets and a grid mesh on Ks such 
that by the performed construction their union is an optimal AM for 



K ■=n (3.5). 



We point out that such a construction seems to be playable also for the 
quite more general class of all bounded connected domains satisfying a UIBC, 
however some not trivial instances should be removed. 

For the reader convenience we recall here the Bernstein Inequality. 

Theorem 3.1 (Bernstein). Let p E ^"(M), then for any a < b E we 

have 

\Pi^)\< n ^ 7=^IHIm- (14) 

\/ [x — a)\h — X) 

Let Q he a bounded connected domain in M'^, we denote by c?[;q(-) the 
distance function w.r.t. the complement of Q and by Tfzni') the metric pro- 
jection onto Cf2. Let us introduce the following notation that can be thought 
as in figure [TJ 

l(x) := min inf Ja > : y + A ~ ^, ^ n\ (15) 
In := inf/(a;). (16) 

We continue to use the same notation as in the previous section for the closure 
and the boundary of Q, namely X := dQ and K := Q. 

The following consequence of Bernstein Inequality will play a central role 
in our construction. 
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Proposition 3.1. Let Q be a bounded connected domain in M'^ and let us 

introduce the sequence of functions 

= , if dznix) < In 



(p^{x) := { V^CnWfe-rfcnC^))' . (17) 

otherwise 



■x-y 
\x-y\ 



For any x E Q let v E {ifzfj '■ V £ 7r[;j^(a;)}, then for any p G ^"(M'^) we 
have 

\\dvP{x)\\ < ipn{x)\\p\\s,{x) < (Pn{x)\\p\\K, (18) 

where we denoted by Sv{x) the segment x + [—d[^Q{x),ln — d[^Q{x)]v. 

If moreover we have /q > 0, let us pick any < 6 < In and define the 
sequence of functions 



ifdzn{x) < 5 



, ^ , otherwise 



Then the above polynomial estimate ([18j) still holds when substituting Lpn,5 to 
(Pn and S^^s ■= x + [-d[^n{x), S - d[^n{x)]v to S^. 

Proof. Pick p e ^"'{W'-). Let us take x e such that d[^nix) < In, then 
we consider the segment Sy{x) where v is as above. The restriction of p to 
this segment is an univariate polynomial of degree not exceeding n, then we 
can use Bernstein Inequality |3.5| to get 



dp{x + ^v) 



n 

< \\P\\S^{X), 



+ dzn{x)){ln - dzn{x) - 



evaluating in ^ = we get 



n 



\dvP{x)\ < =IHl5.(.), (20) 

Vdzn{x){ln - dzn{x)) 

thus the first case of the thesis is proved. 

Otherwise if we take x such that d[^n{x) > In, then B{x, d[^n{x)) C Q and 
hence we can pick a segment of length dzn{x) lying in K and having x as 
medium point. The Bernstein Inequality reads as 

max |9,p(x)|<^|bb(,,,,(.)). (21) 
14 



In particular this is an upper bound for the directional derivative of p 
w.r.t. V. 



Notice also that each right hand former of the above inequalities (20) and 



(21 ) can be bounded using \\p\\k instead of max-norm on the considered balls 
or segments because they are subset of K. 

The last statement follows directly by the special choice of 6 < Iq. The 



r.h.s. former in (19) dominates (case by case) the r.h.s. former in (17) when 



cases are chosen accordingly to (19). 



Actually the above proof proves also the following corollary: it suffices to 
take (19) and substitute , ^ by f in the second case. 



Corollary 3.2.1. Let ^ he an open hounded domain and 6 a positive numher 
such that := {x & VL : d^^^^x) > 5} 7^ 0. Then for any v G S"^"^ we have 

\\dM\K, < -ApWk Vpe ^"(M'^). (22) 



A profitable technique in order to build an AM is to use a norming subset 
of the given compact, thus we introduce the following in the spirit of [19]. 

Let us denote by M'^ the one dimensional Hausdorff measure, that is the 
standard length measure in W^. 

Proposition 3.2. Let Q he a hounded connected domain in M.'^ such that 
Iq > and let < 6 < In. Then 

(i) for any x & Q the map 

is constant, let Fn^s{x) he its value, 
(a) We have 

p (3,) _ |^arcos(l - ifdz^{x) < 6 ^^^^ 

I n ( 7r + In ^^2^1^ j , otherwise. 



In particular Fn^s extends continuously to fl. 
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10 20 30 40 50 

Figure 2: A plot of a section of Fn^s along a segment of metric projection, where S = 10, 
n = l. 

(Hi) Fn^s is constant on any level set of d\^Q{-) and sup^^j^^ F^^s = nn. 

Let us set a^^ := ^ where i = 0,1,... m„ and nin is any positive 
integer greater than 2mT, we denote by F^ ^j the a\^^-level set of F^^s- 

(iv) We have 

K,s = {x e K : dz^{x) = dls} > where 




(v) Let r„,5 := U^ToFj^ 5, then for any p G ^''{R'^) we have 

|b||;,<max{2||p||r„,„|b||A'J. (24) 

Proof. (i) The function ipn,s{-) depends on its argument only by the dis- 
tance function, ipn,s{x) = gn,s{d[;,Q{x)). The length of the segment [y,x] 
is clearly constant when y varies in the set nznix). 

Moreover for any y, z E '^[^q{x) let us denote by Ry^z the euclidean isom- 
etry that maps [y,x] onto [z,x], one trivially has c?[;q(^) = (icn(-Ry,^)(0 
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for any ^ G [y,x\. This is because 'Kzni.i) = V ^'^Y C ^ [^jU] thanks 
to triangular inequahty and thus (icn(0 = \^ ^ vl- 

Thus we have 

ly,x] J[yM 

[y,x\ J[z,x\ 

ipn,5{v)d^^{r])- 

z,x] 

(ii) Let us parametrize the segment a.s y + then we have 

^"•'^''^^ , otherwise. ^^^^ 

The first integral can be solved by substitution: s = |(1 — cos^). The 
integration domain becomes [0,6'^.] where |(1 — cos(6'a,)) = d[^Q{x), while 

the integral itself becomes J^"^ dO = Ox, thus the first of ( [23| ) is proven. 

The second integral is an immediate primitive. Fn,5 depends on x only 
by the distance function, moreover we notice that 

lim arcos fl ^| = 7r = lim f vr + In - 



hence Fn,5 is a continuous function of the distance function, since rfpQ 
is well known to be 1 Lipschitz Fn^s is continuous on f2. 

Since cict^ extends continuously to VL, then Fn^^ do, we must take Fn^slan = 
0. 

(iii) We already used that Fn^i depends on x only by the distance function 
and hence Fn,s\d^ (a) = constant, moreover the functions arcos (l — y) 
and (vr + In l) are both increasing in [0,max^^-^d[^Q{x)], see figure 2| 
hence any level set of Fn^s must coincide with a suitable level set of the 
distance function. 
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(iv) The thesis follows immediately by inverting the equation 



n arcos I 1 — 



2dl 



6 



'n.S 




(v) Let p G ^"(M'^) be fixed, let us pick x E K, then two situation can 
occur. In the first case x & Kg and in this case we have \p{x)\ < \\p\\ks- 
In the second we suppose x ^ Kg, let us consider y G 7i[^fi{x), the 
segment [y,x] cuts \/i : d\g < d^^^x), moreover flT^ ^ = 
thanks to the monotonicity of F„^5 along any segment where is 
monotone. 

Let i{x) := max{i : d\^ < rfcn(^)} ^"^^ ^^t be the (only one) 

intersection of r^^^^"*"^ and the ray starting from x and having direction 
x-y. 

Let s(-) be the arc length parametrization of the segment [|/*(^), |/*(^)+^] 
now we have 



Where in the last two lines we used the special choice of a\^ g as equally 
spaced points in the image of F„ ,5 and the choice of m„ > 2m:. 

To conclude we take the maximum of the above estimates w.r.t. x & K 
thus letting i varying along all 0, 1, ... , m„ — 1 and considering both cases 
X G Ks and x ^ Kg. 

Now we consider the rather regular case when the bounded connected 
domain Q is ^^'^ or, equivalently, preserves a uniform double sided ball con- 
dition. 




71,0 Zi 



< 
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Let us notice that for all ^^'^ domains that has uniform interior ball of 
radius r one has immediately Iq > 2r, we will use twice this trivial conse- 
quence of the regularity assumption without recalling it each time. We stress 
that such estimate does not hold in general for domains satisfying UIBC with 
positive radius: we are fully exploiting the boundary regularity assumption. 

Proposition 3.3. Let Q preserve a uniform double sided ball condition of 

radius r > 6 > 0, let rrin > 2mT, then 

(i) For any z = 1, . . . m„ V]^^ is a ^^'^ hypersurface. 

(a) For any p G ^"(M'^) any x G Tl^ g and any v G S'^^^ HTxTl^ ^ where i = 
0, 1, ... , m„ we have 




\dM^)\<<i:7: ^ (26) 



1,2, . . . ,m„ 



Let 1 < fi and let us pick the finite subsets Y^g C such that 

"'''-^ ' ■ =l,2,...,mn 




(Hi) For any p G ^"(M'^) we have 



lb||r„,. <lbl|y„,, + ^lblk, (28) 



where Y^^s := ^T=oKs- 



Proof. (i) If i = we already known that dQ is a ^^'^ hypersurface from 
the hypothesis of regularity of the domain Q. 

If 2 > it is well known that d[^Q{x) is differentiable at each x G M*^ 
such that 7!'[^q{x) is a singleton [13], thus in all Q \ Ks, at such x we 
have 

Moreover this function is a Lipschitz one when restricted to the set 
Unp{{]n) ■= {x G M'^\C^] : n^nix) is a singleton} D n\Ks [13] [Theorem 
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4.8]. The ^^'^ regularity of the boundary imphes that the Oriented 
Distance function 

bn ■= - dn 

has Lipschitz gradient in a full tubular neighborhood of the boundary 
B{dQ, 6), moreover for any x G f2 \ Ks we have 

x-7ron{x) x~7Tzq{x) 

^bn{x) = — = — = Vd^^ix). 

\x-7ian{x)\ \x-tt[^q{x)\ 

Hence Vcico is a Lipschitz function that extends across the boundary 
of f2 to a Lipschitz function. 

We notice that Vd[^Q{x) ^ implies that Rk V(i[;Q(x) is 1 at any point 
of O \ K^^ therefore any level-set of c?(;q contained in f2 \ ii'^ is a "^^'^ 
d—\ dimensional manifold. 

(ii) Let us pick any x in ^ and let y G '/r[;Q(x), v G T^L^ ^ = Ker(V(i[;Q(x)). 
The ball B^, := B{x + ^''VrfcQ, r') , where r' > 6/2, is a subset of K, 
moreover v G TxB^- 

We recall that on any ball B{xo, p] a Tangential Markov Inequality of 
exponent 1 and constant ^ holds true for polynomials, i.e. for any 
^ G dB{xQ,p\ we have 

\duPm < -|bb(xo,p] g S'^-i n r^dB{xo,p\. 
Therefore it follows that \d^p{x)\ < ^\\p\\b., < ^Ibll^- 



(iii) Now fix any z G {0, 1, . . . ,m„}, by (27) for any a; G F^^ there exist a 
point y E Y^^ and a curve 7 lying in F^ ,5, connecting x to y and such 
that Var[7] < /ipi ^X^n,^) • Let us denote the unitary reparametrization 
of 7 by 7, then we have 

\p{x)\ < \p{y)\+ / -^^m< 



dt 

< Iblly^ + V ,(^n,<5) max < 

".,0 n.A c^T^ . ,,^Sd-ln-y^pi ^ 



< b y», + - b B. < 

„ n 1„ „ 

< \\p\\y^, + -\\p\\k- 
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Let us take the maximum w.r.t. x varying in F^^j and i varying over 
{0, 1, . . . ,m„}, we obtain ||p||r„,, < |b||y„,, + j;\\p\\k- 

It is worth to stress that the combination of previous proposition actually 
gives the following version of Bernstein Multivariate Inequality. 

Theorem 3.5 (Bernstein Multivariate Inequality). LetQ be a bounded 
connected domain in M.'^ satisfying the UIBC with parameter r > 0. Let 
< 6 <2r, n eN and p E ^"'{W^), then we have 



\dvP{x)\ < 



n 



K^v e S'^'^ where 



otherwise. 



Proof. Let x ^ K \ Ks. Let us consider the orthonormal basis 
{h',rii,ri2, . . . jTjd-i} such that u := Vd^^Q^x). Then one has 

\dMx)\ < \{v,Vp{x))\ < |w|i|Vp(x)U < v^|^^|2|Vp(x)U < v^|Vp(x)|, 

where both |Vp(x)|oo and |f |oo are calculated w.r.t. the basis 



(29) 



Now by (\18v and (26) we have 



\duP{x)\ < 



n 



Vdcn{x){S - d^^ix)) 
2n. 



\P\\k 



therefore 



\Vp{x) 



< 



n 



^Td 



\\y{^ di^[x)[b - di^[x)), 5/2} 



mm 



Thus by ( 29 ) we have 



\dvP{x)\ < 



n 



^HVdcn{x)iS - d^^{x)), 6/2} 



\p\\K\/xeK\Ks 



mm 



If on the other hand we pick x E Ks then (|22j) implies directly the second 
case of the thesis. 
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We are now ready for the main result of this section. As for Theorem 2.3 
the proof is fully constructive. 



Theorem 3.7. Let Q be a bounded ^^'^ domain in M.'^, then there exists an 
optimal admissible mesh for K := Q. 

Proof. Since f2 is a bounded ^^'^ domain it satisfies the uniform double 
sided ball condition of radius r > for a suitable r. We fix 6 < r. 

First we build an AM for K, then we will show optimality, let p G ^"(M'^) 
be fixed. 

Let us consider for any A > 1 a mesh An^s,x '■= Yn^s U Zn,5,x such that 

The finite set Zn,&,\ C Ks can be pickecj^ using any uniform grid having step 
size not exceeding where w is the quasiconvexity parameter of K^. That 
is oo > w > 1 such that for any Xi, X2 € there exist a rectifiable arc 7x1,3:2 
lying in Ks and connecting xi to X2 such that Var[7j.^ j^j] < w\xi — X2I 

The existence of such a constant is guaranteed by the regularity of dK^ := 
that is a '^^'^ hypersurface. 

By the grid construction we can bound the cardinality of the mesh by 
diameter: 

Ca.d(Z„«)<(!'!^^)' = 0(„<'). 

\ wXri J 

Now recall that any finite length curve 7 can be re-parametrized by arc- 
length on the interval [0, Var[7]] obtaining a 1 Lipschitz parametrization, thus 



"'Equivalently one can build Zn^s,\ such that the conclusion (30) holds true and having 
the right cardinality asymptotic performing the same construction as in (^[Th.5]. For, let 
us notice that since is smoothly bounded it preserves a Markov inequality of exponent 
2. 
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in particular we can always suppose wlog |7'| =jt^-a.e. 1- By (22) we have 



< max < + min < / 

< \\p\\z„sx + hKs{Zn,5,x) max \d^p{0\ < 

< \\P\\Z„,S.X + hKs{Zn,5,x)-^\\p\\K■ 

\\p\\K,<\\p\\z„,,x+l\\p\\K. (30) 

Since the geodesic fill-distance assumption made on Y^^s satisfies the 
step-size bound required to apply (28), then 

lblk.<lbl|y„, + ^lblk, (31) 

here we must consider /i > 2. 

By the special choice of 6 < r < 1^/2 we can use (24) jointly with the 
above estimates and we obtain 

IblU < max{2||p||y„_, + 2^|b||^ , ||p|U„,, + \\\p\\k}. 
By the elementary properties of max we have 

WpWk < max{-^ , j^}\\p\\y„,suz„,s,x- (32) 

Thus Yn^s U Zn^s^x =: An^s satisfies 

\k < C(5, A,/i)||p|U„,, \fp e ^"(M") \fn e N. (33) 



We are left to prove the optimality of An,s, since we already noticed that 
Caid{Zn,s,\) = 0{n^) we are reduced to prove Card(F„,5) = 0{n^). 

When i = Proposition 2.2 assures the existence of such an Y^^ with 
hvl.iyls) < i and Card(y„°5) = 0{n'^-^). Let us study the case i > 0. 

Let us denote by i/(-) is the outward normal unit w.r.t. F^^ defined as 
follows. First one define u : d*Q — ?■ as the measure theoretic outward 
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normal unit on the reduced boundary d*Q of Q (for these Geometric Measure 
Theory definitions we refer to [IB]), then one can extend z> to a Lipschitz 
function u : dfl — )■ having Lipschitz constant < L < cxd. The definition 
of d*fl takes sense because is a ^^'^ domain, hence in particular a set of 
locally finite perimeter or a Caccioppoli Set, moreover we have dQ \d*fl = 
by boundary regularity. On the other hand the Lipschitz extension of z> takes 
place thanks to [1] [Theorem 1.1]. We know by the Structure Theorem for sets 
of locally finite perimeter ([HI Section 16] or [lOj) that dfl is the countable 
union of compact pieces of manifold Sj and we have j^ls ix) G (TxSj)'^ 
for any x G dQ and a suitable choice of j = j{x). 

Let us notice that is globally defining dQ, is differentiable in a neighbor- 
hood of and Rk Vbn{x) = 1 for any x G d^, thus we have V&n(a^) £ (TxSj)'^ 
and since both vectors are unitary VbQ{x) = ±z/(a;). Finally we notice 
that for any x G dfl one has Vfen(x) = \imQ^y^x^d[^n{y)^ hence — z/ is 
the continuous extension of VcicQ to dQ and trivially x — dl^^u^x) G F^^, 

Let us introduce the family of maps /, := ^vrc^lpi ^ j 

f . -pO , -pi 

Ji ■ ^ n,5 ^ n,5 

X I y X d,^^y{x^. 
{fi}i=i^2,...,mrr is an equally continuous family of Lipschitz constant 
, max {l + Ld{Tl,Xn.6))<{^ + L6)- 

i=l,2,...,m„ 

Now Area Formula says that /j (being 1 + L5 Lipschitz) maps a mesh 
of F° 5 with geodesic fill distance onto a mesh in F^ ,5 having geodesic 
fill distance bounded by h. We already used this property and explained its 



application with more deta ils i n the proof of Theorem 2.3 



Thanks to Proposition 2.2 we can pick the mesh Y^^ C F^^ such that 
V,,(tV) < 2lld^) cardinality bound Card(y;V) = 0{{lf-') 

2M1T5L) by let us set Y^^ 



where we denote „ (^^+xt\ by h, let us set Y^^ := {fi{y),y G Yn^s}, we are left 



to prove 



rrin mn / d—l\ 

Card(F.,,) = J2 Cardr;, = n^'^ + E ^ ' = ^(^')- 

i=0 i=l ^ ^ 

Since Yn^s meets both the cardinality and geodesic fill distance require- 
ments we are done. 
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